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I. Summary 


This project focused on the generation and suppression of oscillatory thermocapillary 
convection in a thin liquid layer. While the bulk of the research was experimental in nature, 
some theoretical work was also done. 

The first phase of this research generated, for the first time, the hydrothermal-wave 
instability predicted by Smith and Davis in 1983. In addition, the behavior of the fluid 
layer under a number of conditions was investigated and catalogued. A transition map for 
the instability of buoyancy-thermocapillary convection was prepared which presented 
results in terms of apparatus-dependent and apparatus-independent parameters, for ease of 
comparison with theoretical results. 

The second phase of this research demonstrated the suppression of these hydrothermal 
waves through an active, feed-forward control strategy employing a C0 2 laser to selectively 
heat lines of negative disturbance temperature on the free surface of the liquid layer. An 
initial attempt at this control was only partially successful, employing a thermocouple 
inserted slightly below the free surface of the liquid to generate the control scheme. 
Subsequent efforts, however, were completely successful in suppressing oscillations in a 
portion of the layer by utilizing data from an infrared image of the free surface to compute 
hydrothermal-wave phase speeds and, using these, to tailor the control scheme to each 
passing wave. 

One Ph.D. student (R. Jeffrey Riley) completed his degree at Georgia Tech under support 
from this grant. Another student who had been working on the project with additional 
support from a NASA GSRP fellowship had to leave the program early due to failure to 
successfully complete the Ph.D. qualifying examinations. Two additional doctoral students 
(Simon Benz and Peter Hintz) from the laboratory of Professor Dietrich Schwabe at 
Universitat GieBen each spent six months in Atlanta working on this project Their 
research on this project will become part of their doctoral theses at Universitat GieBen and 
Professor Neitzel is a member of their doctoral committees there. 

Two papers have been accepted for publication in the Journal of Fluid Mechanics and will 
appear in 1998. Corrected preprints of these papers are included as an appendix. The final 
paper will be a discussion of the results of our experiments in light of several linear- 
stability analyses (including our own) of the basic state of steady buoyancy-thermocapillary 
convection. This paper will be submitted for publication in calendar year 1998. 
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Instability of thermocapillary-buoyancy 
convection in shallow layers. Part 1. 
Characterization of steady and oscillatory 

instabilities 

By R. J. RILE Yf AND G. P. NEITZEL 

The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology 

Atlanta, GA 30332-0405, USA 



(Received 5 October 1996 and in revised form 12 November 1997) 

Combined thermocapillary-buoyancy convection in a thin rectangular geometry is 
investigated experimentally, with an emphasis on the generation of hydrothermal- 
wave instabilities. For sufficiently thin layers, pure hydrothermal waves are observed, 
and are found to be oblique as predicted by a previous linear-stability analysis (Smith 
A: Davis 1983). For thicker layers, both a steady multicell state and an oscillatory 
state are found to exist, but the latter is not in the form of a pure hydrothermal wave 


1. Introduction 


Thermocapillary convection is bulk fluid motion driven by temperature-induced 
mterfacial-tension variations. Such flows occur in situations with liquid-liquid or 
liquid-gas interfaces which also experience variations in temperature. The case of 
interest is that of a free surface between a liquid and a passive gas. For liquids of 
interest here, surface tension decreases with increasine temperature, and this effect 
is represented by the fluid property >■ - -dc/dT, where a is the surface tension 
and 7 is the temperature of the interface. Through the tangential-stress balance at 
the interface, surface-temperature gradients generate an interfacial shear stress which 
drives the surface flow in the direction opposite to that of the surface-temperature 
gradient. For thermally driven flows, buoyancy forces will also be present, but if 
the dimensions of the flow are small enough or if the flow occurs in a microgravity 
environment, thermocapillary forces will be dominant. Examples of such flows include 
the tnermocapillary migration of small bubbles or droplets in a thermal field and 
classical Rayleigh-Bcnard convection. 


Recently, there has been interest in the role of thermocapillary convection in 
processes such as thin-film coating and crystal growth from the melt. Specifically 
“ crystal-growth flows there is interest in the instability of steady thermocapillary’ 
convection, leading to oscillatory convection which results in poor crystal quality 
/rTJ ght , b >, unst ' ad y solidification at the liquid /crystal (freezing) interface 
(Gatos 198.), To this end there has been a large amount of research done on a model 
of tne float-zone crystal-growth process known as the half-zone. In this model a liquid 
bridge is held in place by surface-tension forces between two coaxial cylindrical rods 
maintained at different temperatures. Buoyancy is minimized by heating the half-zone 

+ Present address: Cabot Corporation, Research and Development, Pampa, TX 79066, USA. 
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Figuri 1 . Problem geometry and coordinate specification. 


at the top and cooling at the bottom to provide a stable axial-buoyancy gradient 

lot ° f j he Lq u d u ndge iS ‘ huS SUbjeCl t0 3 tem P erature gradient and the 

he bulk is driven by the resulting thermocapillary convection. 

The first documentation of oscillatory thermocapillary convection in half-zones is 
due to Schwabe er a!. (1978) and Chun & Wuest (1979) Both groups measured the 
oscillating temperature fields in h^Hnel witETiaiMters of a few millimetres and 
found oscillations with frequencies of the order of 1 Hz. More in-depth experimental 
studies of oscillatory flow in half-zones can be found in the work of Preisser Schwabe 
& Scharmann ( 983) and Velten, Schwabe & Scharmann (1991). Recent theoretical 
^ lhe problem has been performed by Shen et al . (1990) Neitzel 

e ‘ ai. (1991, 199^). Uanschura et al (1995) and Levenstam & Amberg (1995) The 
first two of these studies employ energy-stability theory, the next pair, linear-stability 
theory, and the final one, direct numerical simulation of a model half-zone problem 
The results of the various theoretical approaches are in agreement with each other 
and in reasonable agreement with available experimental results (see e.g. Velten et al 
1991)) on half-zones. 

The flow examined in this work is that of a thin liquid layer of depth d in a 
horizontal rectangular slot driven by both thermocapillary and buoyancy forces as 
shown in figure 1. This flow possesses some of the central features of the half zone 
such as a return-flow basic state, a velocity profile for which is sketched in the figure 
The geometry allows for very large aspect ratios L/d and W/d to minimize the effects 
of endwalls and sidewalls, respectively. By changing the depth of the layer, control 
o\er the influence of buoyancy is also possible. A temperature gradient is produced 
on the tree surface by maintaining a temperature difference AT = T H - T > 0 
between the two endwalls; the sidewalls and bottom are assumed adiabatic. Since 
' n n'" e thermocapillary stress pulls the free surface from the hot wall to the cold 
wall Buoyancy also serves to drive the flow in the same direction as thermocapillary 
lorces. but the case of greatest interest in this study is that for which thermocapillary 
enects are dominant, implying thin liquid layers. 

The strength of thermocapillary convection is characterized by the Maraneoni 
number, £ 



Mai ~ 


yd^AT 

fiotL ’ 


(1) 


where n is the dynamic viscosity, a is the thermal diffusivity, and L is the streamwise 
domain length indicated in figure 1. This definition facilitates the comparison of 
experimental results from different length slots, since the driving force for the flow 
ts AT /L; however, it differs from that used in theoretical treatments for layers of 
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infinite horizontal extent, i.e. 

A/ „ yd'dT/dx 

Ma = — ’ ( 2 ) 

which employs the temperature gradient along the interface. We shall refer to the 
quantity defined in (1) as the laboratory Marangoni number, and indicate it by the 

u 131 iP W w g l ° w C existcnce of thermal boundary layers near both the hot 
and cold walls, Ma L > M a where the value ofdT/dx employed in (2) is that in the 

hneTr Th C e°dtfinV 0n ’ °/x^ flow - w . here the surface-temperature variation is nearly 
T he r»t>on of Ma provided in (2) permits easier comparison with available 
r> or infinite layers, but requires an a posteriori measurement of dT/dx. If one 

the ,n . st f a . blht >’ °J tbjs core flow is responsible for observed escalations in 
wry thin layers and that end effects serve merely to modify these, then Ma is the true 
Marangoni number of interest since it is apparatus-independent. As layer depth d 
decreases for fixed I (see figure 1), the core region of the flow exhibits a more nearly 
linear variation of surface temperature with stream wise distance. Consequently, as 
d/L — 0, we expect better agreement between experiment and theory due to the 
existence of a more linear thermal field and also to the reduced effect of the endwalls 
on the core velocity field. 

Buoyancy effects can be quantified in terms of the dynamic Bond number, 
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(3) 


L^hich is a measure of the relative strength of buoyancy forces to thermocapillan 
torces Here p is the density, g is the gravitational acceleration, and 0 is the thermaf- 
expansion coefficient The dynamic Bond number can be formed as Ra/Ma , where Ra 
is a Rayleigh number which uses the same horizontal A T as the Marangoni number. 

Also of importance in these flows is the Prandtl number, Pr = v/a, where v is the 
kinematic viscosity. The Prandtl number for the liquid used in this study is 0(10). 

implies a strong coupling between the thermal and momentum fields. Either of 
the .larangom numbers defined above may be expressed as the product of the Prandtl 
number and a corresponding Reynolds number, i.e. Ma L = Re L Pr or Ma = RePr 
Since the thermocapillary-driven surface flow in the closed domain of figure 1 
creates a pressure gradient to drive the return flow, the normal-stress boundary 

’u 11 ° n the free surface rec l uire s a non-zero surface curvature (Davis & Homsy 
19S0J. However, in the limit of infinite mean surface tension the surface becomes 
flat, and hence a constant-depth thermocapillary slot flow as sketched in figure 1 is 
possible Two dimensionless groups that indicate the amount of surface deflection are 
the surface-tension number 

pdo 

5 ■ -tf’ (4) 


and the capillary number, 


Ca = — . 


(5) 


A flat free surface is achieved if S - cc, and Ca 0. In the present experiments, 
S and Ca are typically of the order of 25000 and 0.001, respectively, indicating that 
the assumption of a flat free surface is reasonable. 

Pure thermocapillary convection and its counterpart including buoyancy have 
received attention from both theoreticians and experimentalists. Sen & Davis (1982) 
computed states of steady two-dimensional thermocapillary flow in slots, including 
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free surface deformation. Smith & Davis (1983) investigated the stability of a return- 
flow basic state in the core of such a slot to convective instabilities, finding a new form 
of instability which they termed a hydrothermal wave. This mode is a temperature- 
disturbance wave that propagates in a direction which depends on the liquid’s Prandtl 
number. For small Pr, the waves propagate transversely across the layer and for 
large Pr the direction is nearly opposite that of the surface flow. For intermediate 
Pr, the hydrothermal waves arc predicted to propagate obliquely against the surface 
flow, and this behaviour is the expected one for the experiments presented in this 
paper. The physics of the instability mechanism were later elucidated by Smith (1986) 
for the small- and large-Pr limits in terms of the interaction between the thermal 
field and the inertially (low Pr) or viscously (high Pr) dominated flow field. For 
the Prandtl number (approximately 14) of the liquid used in these experiments, the 
high-Pr mechanism is the more relevant one. In this limit the hydrothermal-wave 
instability is sustained by the coupling between the convection of thermal disturbances 
b) the basic-state velocity field, and conduction of these disturbances to create new’ 
instability-initiation spots. Extension of the Smith & Davis theoiy to account for the 
influence of buoyancy forces has been done by Garr-Peters (1992a, b), Parmentier, 

Regnier & Lebon ( 1993 ) and Mercier & Normand ( 1996 ). 

Kirdyashkin (1984) was the first to analytically investigate combined thermo- 
capillary-buoyancy-driven flow in a thin two-dimensional slot in the core region 
away from the end walls. He presents both velocity and temperature profiles, and also 
includes experimental results obtained using ethyl alcohol in a 900 mm slot with 
a typical layer depth of 10-20 mm. Steady-state profiles are fully developed in the 
core region, agreeing well with analytical solutions. A parallel-flow core region is 
typically observed in the middle 90% of the domain, with endwall influence extending 
approximately three layer depths from each wall. 

\illers &. Platten (1992) conducted experiments and performed two-dimensional 
numerical simulations of combined thermocapillary-buoyancy flow in a thin layer of 
acetone (Pr = 4.2). A centra! finding of their work is that, as the Marangoni number 
increases, the flow transitions from a single steady convection cell, to multiple steady 
cells, and then to an oscillatory state with oscillation periods on the order of 5 s. De 
Saedeleer ei al (1996)] have performed more recent experiments in decane (Pr = 15), 
observing also that the first transition is to a steady multicellular state. We shall see 
from the present results that a combination of buoyancy and geometric factors are 
likely to be responsible for this deviation from the Smith & Davis (1983) transition 
scenario. 

Schwabe et a!. (1992) have documented similar behaviour in slot and annular- 
pan geometries employing layer depths of ethanol (Pr = 17) on the order of a few 
millimetres Here, multiple steady multicellular structures are also observed. Transition 
Marangoni numbers for the various regimes are determined; that from a single steady 
cell to multiple steady cells occurs for Ma > 600, and the transition from multiple 
steady cells to an oscillatory state is observed for Ma > 3000. Most of the reported 
results are for the annular-pan geometry, for which no theoretical results are presently 
available Schwabe et al. (1992) report two different types of oscillatory flow states, 
depending on the layer depth: for d < 1.4 mm, a short-wavelength azimuthally 
travelling state is seen, while, for d > 1.4 mm, a long-wavelength state with a radially- 
propagating component is observed. Results of the present work could imply that the 
difference between these two states may be related to an effect of the dynamic Bond 
number. Unfortunately, very few- results are reported for the rectangular geometry, 
which is very similar to that employed for the present experiments. Those which are 
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provided are for the single layer depth of 1.8 mm, which is likely to be too deep to 
observe pure hydrothermal waves, as will be noted later. 

"^ e Primary aim of the experiments described in this work is to examine com- 
bined therm ocapillary-buoyancy convection in thin horizontal layers, with the goal 
of observing a pure hydrothermal-wave instability. This has been successfully ac- 
complished. Since this instability is believed to be responsible for the degradation 
of material grown via the float-zone process, the possibility of eliminating such os- 
cillations is of interest. In Part 2 (Benz et al. 1998), the feasibility of suppressing 
hydrothermal waves through active periodic heating of lines on the free surface is 
demonstrated for the moderate-? r liquid used in these experiments. However, it must 
be emphasized that the hydrothermal-wave instabilities which occur in a liquid of 
moderate Prandtl number are driven by a different mechanism (Smith 1986) from 

that responsible for similar instabilities in low-Prandtl-number liquids characteristic 
ol liquid metals. - 


2. Experimental setup and techniques 

2.1. Working liquid and apparatus 

The working fluid used in these experiments is 1 cS Dow Corning silicone oil. It 
has the desirable properties that it is transparent to visible light, allowing both flow 
visualization and quantitative measurement using laser-DoppIer velocimetry (LDV), 
and has a free surface which is relatively resistant to contamination. In addition, we 
shall examine the spanwise structure of the observed instabilities by visualizing' the 
thermal field of the free surface using an infrared camera. The chosen oil also has 
an absorption band in the infrared within the detectable wavelength range of this 
camera which is strong enough to allow infrared visualization of the liquid’s surface 
temperature rather than that in the bulk. 

The Prandtl number of the silicone oil used here is Pr = 13.9, indicating that the 
thermal field is strongly coupled to the velocity field and, according to the theory 
of Smith & Davis (1983), the preferred instability should be a hydrothermal wave. 
The properties of the oil, coupled with the thin layers investigated here (0(1 mm)) 
also ensure that the free surface of the layer is flat. A fluid possessing a low vapour 
pressure is also desirable in thermocapillary flow experiments due to the likelihood 
of evaporation at the hot end of the apparatus. The 1 cS oil was chosen after 
experimenting with silicone oils with viscosities in the range of 0.65-10 cS in part 
because it was found to have the fewest problems with evaporation. 

The test cell w as designed with the goal of observing a hydrothermal-wave insta- 
bility: a schematic of the cell is shown in figure 2. The length L and width W of 
the flow section are 30 mm and 50 mm, respectively; for the range of layer depths 
examined in this work, these choices afford streamwise (x-direction) and spanwise 
(y-direction) aspect ratios of 12 < L/d ^ 40 and 20 < W/d 67, respectively. These 
permit both the establishment of a well-developed return-flow velocity profile in the 
core region aw ay from the walls and the accommodation of a relatively large number 
°[ hydrothermal waves. For a liquid w>ith Pr = 13.9, the linear-stability analysis 
of Smith &. Davis (1983) predicts a most-dangerous disturbance with a wavelength 
of approximately 2.4 d, propagating upstream at an angle of 20 c to the streamwise 
direction of the basic state. Thus, at least a dozen hydrothermal waves can exist in the 
domain, and the large spanwise extent of the apparatus allows minimized impedance 
in the oblique propagation direction. 
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(b) 


Machined lip 



Experimental apparatus schematic (all dimensions are in mm) 
£{b) side view. 




top view; 


To provide an accurate determination of the depth d and its continuous variability, 
the apparatus was designed with a movable base and reservoir, as seen in the side view 
of figure 2(b). The base is a machined Plexiglas block which fits snugly within the cavity 
created by the endwalls and sidewalls, yet has very small gaps (0(0.025 mm)) between 
it and the vertical walls which permit the addition of liquid to the layer through 
the reservoir Three ball bearings mounted in the underside of the base ride on 
micrometers, allowing depth adjustment of the flow domain to within 0.005 mm. The 
apparatus is filled by injecting sibcone oil into the reservoir, as suggested above. The 
sidewalls are also constructed of Plexiglas in order to approximate adiabatic boundary 
conditions The endwalls, through which temperature-controlled water is circulated, 
are aluminium to provide uniform constant-temperature boundary conditions at both 
endwalls Endwall temperatures are measured directly using thermocouples mounted 
in the centre of each endwall a distance of 1 mm from the surface which is in contact 
with the oil The difference in endwall temperatures is able to be resolved to within 
0.05 C. 

In order to achieve a definite and repeatable layer depth without a meniscus, a 90 c 
2 mm lip has been machined along the top of the sidewalls and endwalls which serves' 
to pin the contact bne of the silicone oil. The upper surface of this bp is coated with 
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a fluorad (3M Scotchguard) surface-modifying treatment, which inhibits the liquid 
from wetting the lip. The existence of this bp with a pinned contact line is crucial to 
the establishment of a layer of uniform depth. Levelbng of the apparatus is achieved 
bj mounting it on a three-point micrometer-type tilt table. Adjustments to the tilt are 
made while observing the reflections of a grid in the free surface near the lip regions, 
and continue until the grid appears undistorted at all points along the lip, indicating 
a flat free surface. Isolation from laboratory vibrations is achieved by securing the 
tilt table to a commercial vibration-isolation table. 

Two approaches were investigated for minimizing disturbances at the free surface 
due to transient laboratory air currents. The first was to cover the apparatus with a 
Plexiglas lid allowing a 5 mm thick gas layer above the free surface. However, this 
method experienced problems with condensation of trapped sibcone-oil vapour at the 
cold wall The formation of droplets on the cold horizontal lip ultimately results in the 
loss of the pinned contact bne, requiring the cessation of experiments and the cleaning 
and re-coating of the apparatus. An alternative technique was successful, however. A 
Plexiglas shroud 60 mm x 60 mm in cross-section and 100 mm in height was placed 
around the flow cell. Since the upper end of the shroud was open to the laboratory, 
over time liquid depletion due to evaporation results in a change in the position of the 
free surface. Although the rate of evaporation is quite small (0(0.1 pl/s)), the resulting 
depth change is unacceptable for these experiments. Hence, the evaporating sibcone 
oil was replaced by liquid supplied from a constant-head container to the reservoir 
at the required rate. At the typical re-supply rate, the effective vertical velocity at the 
free surface is roughly 0.07 pm s * or about 0.001% of the measured free-surface 
speed, and is therefore felt to have a negligible impact. 

Properties of the Dow Coming 200, 1 cS silicone oil used in these experiments were 
obtained from the manufacturer, with the exception of the surface tension. These 
data were measured directly by colleagues at the Microgravity Advanced Research 
Support (MARS) Center in Naples, Italy using the wire loop method for the range of 
the temperatures encountered in the present experiments. From these measurements, 
a linear least-squares fit of the data yielded the result 

<r( T) = 1 7.237 - 0.0755( T - T 0 ) mN m"' (6) 

where T 0 = 25 ; C was the reference temperature about which experiments were 
performed. 


2.2. Experimental techniques 

Four different experimental techniques were used to explore the flow’ states encoun- 
tered in this work. Quantitative measurements of steady-state velocity fields, both 
prior to and subsequent to the onset of instability, were performed with a TSI single- 
component LDV The flow was seeded with dilute amounts of sibcon-carbide particles 
with a mean diameter of 1.5 pm. Velocity-profile measurements through the depth (z) 
of the layer allow one to see the degree to which the assumed return-flow- basic state 
is achieved; measurements as a function of streamwise (x) and spanwise (y) distance 
permit an assessment of the influences of endwalls and sidewalls, respectively. 

Flow. visualization experiments were performed using two techniques. For the first, 
the floff was seeded with polystyrene microspheres with diameters ranging from 1- 
1 5 pm and these were observed within a 1 mm thick laser-sheet illumination of an 
(x, z)-plane A Dage-MTI CCD camera was used to record the resulting images 
for both steady and unsteady flow conditions. Since the particles were slightly more 
dense than the test bquid, problems were experienced with keeping the particles from 
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a stat ic fluid, particle settling time through the depth of 
the layer (at 10 pm s for an average 8 pm particle) is roughly two minutes; for 

for" thl ,C m^T n , ent K part r c ! es r remain “ sus P ension longer, but insufficiently long 
[* h th °d to be useful for extended periods of time. The second technique 
employed was shadowgraphic visualization across (i.e. in the y-direction) the layer 
Thermally induced mdex-of-refraction variations, integrated through the span of the 
ayer were likewise recorded using the CCD camera. While no attempt was made 
to extract quantitative data from these results, the shadowgraphic technique is quite 
sensitive, providing clear indications of both the existence of multicellular structures 
and unsteadiness. It is the visualization method of choice for pinpointing the flow 
transitions to be reported in §3. The only limitation of the technique, with regard to 
i t, P ? experiments, is the fact that the mean temperature gradient bends the 
hght downward on its way to the camera so that, at some limiting depth, it is no 
longer possible to obtain clean shadowgraphic images of the flow. For this apparatus 

ail Ti. the c COndltIOnS ° f theSe ex P eriments - this limiting depth is 0.75 mm. 

The final technique used in these experiments was infrared thermography of the 
free surface. The use of this technique was central to the investigation of oblique 
hydrothermal-wave instabilities and to the suppression of these waves reported in 
Fart r:. lTi add,tl0n > 3 * * * * * * * 11 permits the observation of the structure of the multicellular 
instability’ and its subsequent transition to oscillatory flow. These measurements were 
performed using an Amber AE-4128 infrared camera employing a 128 x 128 element 
lndium-antimonide focal-plane-array detector, which is sensitive to radiation in the 
range of 1-5.5 urn. The output from the camera was captured with a variable-scan 
frame grabber mounted in an IBM-compatible 386 computer, allowing it to be 
displayed on a video monitor or recorded onto video tape. Further details on this and 
all other techniques employed in these experiments may be found in Riley (1995). 


3. Results 

The combination of methods described above was used to investigate the flow 

structures observable in the test cell over a range of conditions which was as large 

as possible, given the inevitable constraints imposed by the apparatus, test liquid 

and measurement techniques, some of which are described above. The results to 

be p reseme d here will focus on the return-flow basic state, its instability to both 
hydrothermal waves and to a steady multicellular state, and the onset* of time- 
dependent flow from the steady multicellular state. The transition map for the states 
observed is presented in terms of the laboratory Marangoni number Ma L (equation 
(1)), and true Marangoni number Ma (equation (2)) in figure 3. 

There are three major categories into which the observations fall: steady unicellular 
flow (SUF); steady multicellular flow (SMC); and oscillatory flow. The last of these 
categories, however, may be further divided depending on whether the state preceding 
the onset of oscillatory flow was unicellular (leading to the hydrothermal waves (HTW) 
predicted by Smith Sz Davis 1983) or multicellular (leading to oscillating multicellular 
flow (OMC)). It should be pointed out that, while we have chosen to characterize the 
transitional Marangoni numbers in terms of the dynamic Bond number there is the 
additional possible dependence of the transition sequence on the two aspect ratios 
L/d and W/d. Since this apparatus does not allow these to be held constant as the 
depth of the layer is varied, it is not possible to unequivocally attribute the observed 

transition to multicellular flow to a particular effect. 

As noted in §1, Ma L exceeds Ma, owing to the existence of thermal boundary 
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NorTnaIlz ' d IR ca ™ ra emission for the steady unicellular state (just prior to onset) alone 
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ayers at the temperature-controlled endwalls, particularly for moderate-to-high-Pr 
liquids such as the one employed here. These boundary-layer effects become more 
pronounced as the layer depth increases, as seen in figure 4, which shows a normalized 
emission from the IR camera along a line in the middle of the span of the laver as 
a function of streamwise distance x for two different layer depths. For both cases in 
ngure 4, the emission profiles correspond to a state of steady unicellular flow just prior 

° f CIther hydr0thCrmaI WaVCS (figUrc 4a) or stcad y multicellular 
In order to correct the Marangoni number for this effect, the emission data in 
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^ n i ,'° n 1 num , bers versus depth and dynamic Bond number for transitions 

Sterielr?™? waves (HTW), steady multicells (SMC) and oscillating multicells (OMC). The 
mlest gatTd * ““ W “ n0t 3 true transition 10 HTW, n.o. denotes not observed, and n.i. not 


fhe interior region of the flow, i.e. 5 s$ x < 25 mm, are fit with a straight line, 
which, along with the calibration data (including the known measured temperatures 
of the endwalls). is used to extract dT/dx for the core region. This same procedure is 
employed for determining the core temperature gradient used to correct the transition 
map for the steady -» oscillatory' multicellular regime as well, although the data in 
this case are of the type seen in figurd^fa ). In all cases, the temperature gradient is 
determined for a state just prior to the onset of the transitional state. De Saedeleer 
eI al ( 1996 ) also recognized the importance of determining the core temperature 
gradient, which they did using a thermocouple inserted below" the free surface. 

_ One advantage to the present method for determining the true temperature gradient 
C J /c * occurring in the core of the layer is that there is no need to determine 
the effective emissivity of the liquid or to calibrate the grey-level values against a 
blackbody emission source By taking the ratio of the linear fit to the grey-level data 
in the core to the gradient represented by the maximum and minimum grey levels 
occurring at the endwalls of the same data, the absolute values of the free-surface 
temperature need not be determined. 

Since the determination of Ma requires an i posteriori measurement of the core 
surface-temperature distribution, this was not necessarily done for cases not associated 
with transitions between states. Rather, the more convenient laboratory Maraneoni 
number Sla L is reported in the following sections when comparing to specific" ex- 
perimental cases. The transitional values of both Marangoni numbers are reported 
as a function of both the dimensional depth and dimensionless Bo D in figure 3 and 
table 1. 


3.1. Steady unicellular flow 

Of interest in characterizing this state is determining the degree to which the apparatus 
utilized for these experiments permits the establishment, in the core region away from 
vertical boundaries, of the return-flow basic state assumed in the theoretical analysis 
of Smith & Davis (1983). To this end, detailed LDV measurements have been 
performed for two states. The first case, with ( Ma L , Bo D ) = (350, 0.142) corresponds 
to a layer of 1.0 mm depth at a Marangoni number which is yk of t he value at 
which transition to the hydrothermal-wave state was observed. The second case, with 
(Mai, Bop) = (475, 0.320), is for a 1.5 mm layer depth, which first experiences a 
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Figl-R£ 5. Horizontal 'surface' velocity near the hot wall: (a) d = 1.0 mm; (b) d = 1.5 mm. 

transition to steady multicellular flow. The chosen Marangoni number is the same 
fraction j y^/o) of the transition value to the new state as for the first case. 

figures 5 and 6 show the variation of the x-cornponent of ‘surface’ velocity near Ca/^a > 
th*? hot and cold walls, respectively, for the two cases. These measurements were 
performed with the LDV fib^optic head oriented above the layer, looking downward. JtrT 1 
It is estimated (Riley 1995) that the centre of the measurement volume (an ellipsoid j 

with a major axis of 120 pm) is roughly 60 pm below the free surface. Measurements 
for both hot walls show an overshoot of the core velocity followed by an equilibration 
which takes place within a distance x/d = 3 in both cases. At the cold walls no such 
overshoot is observed and the effect of the endwall is constrained to a region of 
approximately x/d = 2/3 in length. 

It is possible that the observed overshoot and recovery at the hot wall are due 
to the formation of a strong roll cell there, in some sense a precursor to the steady 
multicellular state to be discussed in §3.3. Such a cell was observed even for low-Bo 0 
cases at Marangoni numbers immediately preceding the transition to hydrothermal- 
wave instability. The overshoot in velocity is consistent with LDV measurements and 
calculations performed by Villers & Platten (1992) for flows of acetone (Pr = 4.24) 
in containers of modest aspect ratio. 

Another numerical study which has relevance to the endwall region velocity mea- 
surements is that of Ben Hadid & Roux (1990), who performed numerical calculations 
of thermocapillary convection in a long horizontal layer for a low-Prandtl-number 
fluid. Since Pr is assumed to be small, the momentum and thermal fields are de- 
coupled, and the case of a constant applied shear stress at the free surface was 
considered. Thus, there are significant differences between this work and the high-Pr 
situation considered presently. One of the cases computed by Ben Hadid & Roux 
has an aspect ratio L/d - 25, which is comparable to the present 1 mm and 1.5 mm 
layers which have L/d = 30 and 20, respectively. For the lowest-Reynolds-number 
case (Re = 66.7) presented by Ben Hadid & Roux, the free-surface velocity increases 
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x (mm) 


Figure 6 Horizontal ‘surface' velocity near the cold wall; ( a ) d «= 1.0 mm; (b) 4 = 1.5 mm. 


monotomcaliy away from the hot wall and decreases in a similar fashion near the 
cold wall, exhibiting no overshoot at either end of the slot, contrary to what has 
been observed in the present experiments, although this is expected in light of the 
free-surface conditions assumed in the calculations. The two experimental cases under 
discussion correspond to Reynolds numbers of Re = 23.4 and 34.1, respectively, for 
the 1.0 mm and 1.5 mm layer depths. 

Measurements of surface velocity (similar to those of figures 5 and 6) made in the 
vicinity of the adiabatic sidewalls indicate that the free-surface velocity reaches its 
core value much more rapidly than observed near the heated endwalls (Riley 1995). 
For the cases examined, this adjustment at the streamwise centre of the layer takes 
place within a single layer depth. 

In the core region, away from both endwalls and sidewalls, LDV measurements 
made through the layer indicate the degree to which the return-flow basic state is 
attained. The quadratic profile of Smith & Davis (1983) is modified by the existence of 
buoyancy. Solution of the two-dimensional Boussinesq equations for fully developed 
flow in the core yields the horizontal velocity, 

C/( z )»2^_, r + Boc( _i r3 + i z2 _i r)i (7) 

where length has been scaled by the depth d and velocity by the thermocapillary 
velocity 

rr y(dT/dx)d 

v - — i — ■ <*> 

Figure 7 shows measured and theoretical profiles of V(z) for the two layers con- 
sidered in this section. Measurements were performed by having the two intersecting 
beams of the LDV enter the flow field through one of the adiabatic sidewalls. This 
orients the major axis of the ellipsoidal measurement volume normal to the z-direction 
in which the largest velocity gradients occur. Velocity-gradient broadening is there- 
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Figure 7. Basic-state velocity profiles for: (a) d = 10 mm, Ma = 350; open circles, LDV measure- 
ments, solid line, theory for Boo — 0; (b) d = 1.5 mm, Ma = 475; open circles, LDV measurements; 
solid line, theory for Boo = 0.32. 


fore confined to the 40 pm waist dimension of the measurement volume. In order 
to obtain data near the bottom (free surface) of the layer, the plane of the beams 
entering through the sidewall is angled slightly downward (upward) to position the 
measurement volume as close as practicable to the location of interest. Near the 
free surface, where the uncertainty in the measurement is largest, the size of the 
standard deviation in the dimensional velocity for the 1.0 mm layer of figure 1(a) is 
0(0.2 mm s -1 ). compared with a maximum measured speed of roughly 4 mm s _l . 
In the cases shown in figure 7, the velocity has been scaled by its maximum value 
at the free surface, as opposed to the thermocapillary velocity scale. For the thinner 
layer, the comparison between the measurements and theory is made for Bo D — 0, 
indicating the degree to which this layer conforms to the theoretically assumed basic 
state of Smith & Davis (1983). In both cases, the comparison between measured data 
and theoretical profiles is excellent. 

One final question regarding the unicellular state observed in these experiments 
is the degree to which flow in the core region is independent of the x-direction 
of figure 1. In order to assess this, measurements of the maximum value of the 
interior (backflow) velocity were performed at several locations between the two 
heated endwalls, as shown in figure 8. The points displayed as circles are all scaled 
using a constant value of V s (equation (8)) corresponding to the conditions at the 
rightmost cold point, where the properties are computed using a conduction solution 
for the temperature. These data indicate the existence of a fair amount of variation 
between the endwalls. Speculating that a large part of this variation may be due to 
the temperature dependence of viscosity, we have adjusted the velocities depending on 
the local value of U„ computed using the same conduction solution for temperature. 
The adjusted values, shown in figure 8 as squares, show much less variation, lending 
some credence to the hypothesis. 
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Fl °rIfi, M r™ Unl measured backflow velocity versus streamwise distance: circles, scaled by a 
*’ squares, scaled by tA[T(x)], where T is obtained from a conduction solution. 



F' GURf * Instantaneous thermograph of hydrothermal waves on a 1.0 mm deep layer viewed from 

,\ b T:7 h l h0t 3 , nd C0ld Wa!ls are a, the left ' and ngh'-hand sides of the figure, respectively, and 
tne h>arotnermaI waves are propagating from lower-right to upper-left. 


3.2. Hydrothermal waves 

The state of greatest interest in the present experiments is that which results from 
instability of the unicellular basic state at low values of the dynamic Bond number. 
These hydrothermal waves are travelling thermal waves which, for the Prandtl number 
of the liquid used in these experiments, propagate obliquely, with a component which 
opposes the direction of free-surface motion. For Pr = 14 and Bo D = 0, the theory of 
Smith & Davis (1983) predicts a wavelength of approximately 2.5 d and a propagation 
angle ¥ *= 20° with respect to the negative x-direction. In these experiments, the 
hydrothermal-wave transition is observed for Bo D < 0.22, corresponding to layer 
depths d ^ 1.25 mm. 

The onset of time-dependence in this case is an especially sharp one to which 
the shadowgraphic visualization technique is particularly sensitive. The qualitative 
stucture of the flow, when viewed from above with infrared thermography is clearly 
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in agreement with that predicted by Smith & Davis (1983). This is shown in figure 9 
for a la>er of 1.0 mm depth at a condition which is slightly supercritical. In accordance 
with the nonlinear theory of Smith (1988), only one of the two possible families of 
hydrothermal waves is observed. If the apparatus is level and the base and free 
surface are parallel (which is essential to the attainment of this state), either of these 
families may appear. It is clear from this thermograph that the hydrothermal waves 
completely fill the domain with little influence from the endwall boundaries. The same 
is true of the sidewall boundaries, as well, although not observable in this figure; 
they are constructed from Plexiglas, and hence behave as nearly adiabatic boundaries! 
Although very difficult to observe in figure 9, there does exist a single steady roll 
cell located at (and parallel to) the hot wall toward which the hydrothermal waves 
propagate. This cell is consistent with its appearance even in the states which have 
been termed unicellular. It is also apparent from the figure that the wavelength 
increases slightly as the waves propagate across the layer from the cold wall to the 
hot one. This could be a ramification of the fact that the core flow is changing with x 
due to the temperature dependence of viscosity, as discussed earlier, or the fact that 
the wavelength depends on the Prandtl number, which also varies locally with the 
temperature. 

The temperature perturbations associated with these hydrothermal waves may 
be obtained from the infrared camera data, fgure 10 shows both the temperature 
variation (figure 10a) and perturbations (figure 10b) obtained by subtracting the 
mean-temperature variation along a line perpendicular to the wave fronts for data 
extracted from the thermograph of figure 9. From these data, it is seen that the 
amplitude of the temperature perturbation is not constant, but increases in the 
direction of propagation from roughly 0.2 e C to 0.35 °C peak-to-peak. 
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,/t\ble 


d (mm) 

2 (mm) 

/(Hz) 

0.75 

1.95 

0.98 

1. 00 

2.58 

0.57 

0.25 

3.37 

0.38 


c (mm s' 1 ) 

1.91 

1.47 

1.28 


F (deg.) 
27.9 
23.2 
22.4 


2. Dimensional 


wavelength, frequency, phase speed and propagation 
Zjayer depth for hydrothermal-wave cases. 




a function of 



Bop 

0.080 

0.142 

0.222 


Vd 

2.60 

2.58 

2.70 


m 

0.0236 

0.0217 

0.0174 


c/C/, 

0.0614 

0.0561 

0.0469 




TABLf 3. Scaled wavelength, frequency and phase speed as a function of dynamiclBond number 

Jfor hydrothermal-wave cases. 


L 


The structure of the hydrothermal-wave instability is quantified in terms of wave- 
length /. frequency /, wave speed c, and angle of propagation P for the three cases 
examined in this work. All of the quantities were extracted from the infrared images 
of the free-surface temperature fields at values of AT slightly above critical Dimen- 
sional values of these quantities are presented in table 2; dimensionless values of A, / 
and c are given in table 3, where wavelength has been scaled by depth, wave speed 
by the thermocapillary velocity scale U s (equation (8)) and frequency by UJd. The 
dimensionless wavelength is seen to be comparable to the value ).jd = 2.4 predicted 
by the linear-theory zero-gravity analysis of Smith & Davis (1983) and is nearly 
independent of the depth of the layer. The variation of the other quantities with Bo D 
is more substantial, so that comparison with a zero-gravity theory is less appropriate 
A linear-stability analysis incorporating buoyancy has been performed for the Smith 
& Da\is problem (Riley & Neitzel 1998) and a comparison of these results with the 
present experimental ones is discussed in detail there. Suffice it to say here that the 
linear-theory critical Marangoni number, Mac (plotted as the symbol x in figure 3), 
agrees almost exactly with Ma L for the smallest value of Bo D measured in these 
experiments, and that the trend in Ma c with increasing Bo D is the same for both 
theory and experiment. 

^ rev * ous resear chers have measured oscillations in thin thermocapillary-driven 
liquid layers, but have not provided conclusive evidence of the hydrothermal-wave 
instability. The experiments of Villers & Platten (1992) in a thin acetone layer (Pr = 4) 
utilize one-point LDV measurements to determine whether the flow at a single point 
in the layer was steady or oscillatory. In this manner they map out a boundary 
in (Ma, Ra )- space denoting the steady/oscillatory bo u n dary No discussion on the 
structure of the flow producing the oscillations is given, since their measurements 
do not provide that type of information. Much of their data are for relatively thick 
layers (d ^ 2.0 mm, Bo D ^ 0.362), and hence the observed oscillations are probably 
the oscillating multicellular flows to be discussed in §3.4, rather than hydrothermal 
waves. The study by De Saedeleer et al. (1996) reports the first instability in the 
form of a steady structure which then transitions to an oscillatory state. However, 
calculation of the Bond number associated with their experiments reveals that it is 
in excess of three times the maximum value determined for the first transition to a 
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hydrothermal-wave instability in the present experiments. In addition, their spanwise 
aspect ratio was quite small, which may have adverse effects for oblique waves. 

Schwabe et al. (1992) discuss two different types of oscillatory thermocapillary flows 
in a thin annular geometry, so-called short- and long-wavelength instabilities. Their 
short-wavelength instability occurs in layers of d < 1.4 mm ( Bo D < 0.18), and has a 
wavelength of /. as 6 d, for a 1 mm thick layer, whereas / = 2.58 d for the 1 mm layer 
of the present study. The short-wavelength instability is investigated through single- 
point thermocouple measurements, and a shadowgraph technique that illuminates the 
oscillating meniscus formed at the lip of the inner wall of the annulus. The time 
variation of the meniscus at the lip appears wave-like, but no flow visualization in the 
bulk of the fluid is available. The dimensionless wave speed of the short-wavelength 
instability is c = 0.06 for a 1 mm thick layer, which is in agreement with the value of 
c = 0.056 for the 1 mm layer of the present study. Thus, it is possible that Schwabe et 
al. did indeed observe hydrothermal waves in an annular geometry, but the evidence 
is not completely conclusive on this point. 

Ezersky et at. (1993) claim to have produced hydrothermal waves in 5 cS silicone-oil 
layers in the range of 1.2 ^ d < 3.1 mm, but results are presented only for a 3.1 mm 
layer. The flow is investigated using a shadowgraphic technique, and an image is 
presented which is reportedly of a hydrothermal wave. For a layer depth of 3.1 mm, 
Bo d =r 1.4. which is much too large for the generation of pure hydrothermal waves 
according to the present results. Ezersky et al. present a shadowgraph image showing 
a disturbance with a wavelength of /./d » 1.3, which is in the range of the steady 
multicellular flow structure discussed in the next subsection. 

In a similar study. Daviaud & Vince (1993) report an observed travelling-wave 
disturbance in 0.65 cS silicone oil, for layer depths of 0.8 d ^ 2.5 mm. However, 
contrary to the results of the present study, the angle of propagation of the waves 
was found to be perpendicular to the applied temperature gradient (IP = 90' ). This 
could be due to the fact that their rectangular domain is quite large in the spanwise 
direction (200 mm), and relatively short in the streamwise direction (10 mm) in 
which the temperature gradient is applied. They also provide a graph of the critical 
applied-temperature difference (AT„„) as a function of layer depth, which shows 
ATcrr decreasing with increasing depth for 0.8 < d ^ 1.4 mm, followed by an increase 
for depths up to 2.5 mm. This behaviour is inconsistent with the results of the present 
study in which AT m , decreases monotonically (Ma L oc d : ), and could be due to the 
difference in the geometries of the two flow domains. 

3.3. Steady multicellular flow 

The steady multicellular state observed at higher values of the dynamic Bond number 
are spanwise-umform co-rotating rolls. Figure 1 1 shows the structure of this flow as 
observed by the three types of visualization techniques described in §2. In figure 1 1(a), 
a time exposure of sheet-illuminated particles clearly identifies the cellular structure 
and the associated saddle points in a layer of 2.0 mm depth; figure 11(h) shows an 
instantaneous snapshot of a combined shadowgraph/particle image; and figure 11(c) 
provides an overhead view of the temperature field visualized with infrared thermog- 
raphy. This last photograph clearly illustrates the two-dimensional nature of this flow 
state and the relative lack of influence of either sidewalls or endwalls. 

Although it is impossible to ascertain from the thermographic image of figure 11(c), 
the strength of the cells decreases as one moves away from the hot wall. This 
may be observed by examining the streamwise variation of the horizontal velocity 
at a fixed vertical location, as shown in figure 12. This result is for a 1.5 mm 
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(a) 



(b) 



(c) 



S, , eady - mult,ce,lula : flow as visualized by: (a) panicle-path time exposure for a 20 mm 
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Horizontal velocity at z/d ■■ 


■■ 0.9 for steadv/multicellular flow in a 
1^ Mo « 800. 


£ 


(>5 

mm deep layer at 


,ay er depth at a value of Ma L = 800, which is significantly higher than the steady 
multicellular-flow transition value of Ma L - 595 for this Bo D . Note that, even at this 
large temperature difference, the multicells decay rapidly in strength, vanishing at a 
distance of roughly x/L = 2/3. Also of note is the fact that the average velocity 
decreases when moving toward the cold wall, consistent with the earlier observation 
of a temperature-dependent viscosity effect. 

The primary quantitative result to be obtained for this flow state, in addition 
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Figure 1 s. Dimensionless wavelength for steady multicellular stales versus dynamic Bond number. 






Instantaneous shadowgraphs at three instants of time (r is increasin 
bottom) for the oscillatorymuhicellular state. 


r 


from top to 


to the transition Marangoni number, is the dependence of the wavelength of the 
multicellular state on the dynamic Bond number. This is provided in figure 13, and 
exhibits relatively smooth, yet significant ; variation. Once again, it is not possible 
to determine whether these observations are due to buoyancy or aspect-ratio effects. 
However, the degree of smoothness of the results indicates that the observed increase is 
certainly not due to abrupt quantum changes in the number of cells in the apparatus. 


3.4. Oscillating multicellular flow 

If the steady multicellular flow states of §3.3 are driven vigorously enough by increas- 
ing the Marangoni number, a second transition will take place to a f time-dependent 
multicellular state. This state, although oscillatory, is not a pure travelling-wave state. 
This is evidenced by the sequence of shadowgraphs shown in figure 14, which shows 
the flow in a central portion of the apparatus. Time is increasing from top to bottom 
in the sequence, and the hot wall is to the left. Upon examining this sequence, there 
does not appear to be a well-defined wave speed associated with this state; there 
are pronounced differences between what is observed here and similar shadowgraphs 
for the pure hydrothermal-wave state of §3.2 (see Riley 1995). It is difficult to draw 
definite conclusions from such shadowgraph images, however, since the shadowgraph 
is integrated through the spanwise extent of the domain. 

The free-surface thermograph shown in figure 15(a) provides greater insight into 
the structure of this flow. It is clear that there remain strong steady multicellular 
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structures near the hot wall and a pair of oblique waves which appear to propagate 
through them, beginning at the cold wall. This is more evident in figure 15(fc) in 
j IC £ j* 3 a PP ea nng in figure 15(a) have been enhanced through processing with an 
edge-finding filter, and subsequent binary windowing. A chevron pattern is observed 
near the cold wall which is suggestive of two interacting hydrothermal waves. In 
the absence of buoyancy, the nonlinear theory of Smith (1988) states that a single 
family should be observed (such as seen in figure 9) following instability of the pure 
return-flow basic state. However, the oscillatory instability observed here for large 
dynamic Bond numbers is reached from a different basic state (steady multicellular 
flow) than that considered by Smith, so his results are not strictly applicable. The 
image shows that the steady multicells are still established at the hot wall despite the 
vigorous oscillatory convection on the cold-wall side of the domain. 

An oscillatory multicellular state is also observed by Schwabe et al. (1992) in a 
rectangular domain, but the reported results for that apparatus are sparse They note 
however that a 1.8 mm deep layer of ethyl alcohol experiences a transition from a 
unicellular flow to a steady multicellular flow, and subsequently to a time-dependent 
flow at Ma L = 3000. Villers &. Platten (1992) make no attempt to examine the 
structure of their oscillatory flow states in acetone, but for their layer thicknesses 
(d > 2 0 mm, typically), the results of the present study strongly suggest that they are 
actually observing oscillating multicellular flow. 


4. Discussion 

Experiments have been conducted to investigate the hydrothermal-wave instability 
of thermocapillary convection in a thin liquid layer first predicted by Smith & Davis 
(1983). The results of these experiments offer the first conclusive laboratory proof of 
the existence of this instability; results of previous experiments conducted in deeper 
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layers have been clouded by the presence of buoyancy, geometric effects, or both. 
Present results for Boo show a transition to oscillatory flow at a Marangoni 

number of Ma = 324, compared with the theoretical result (for a Prandtl number 
corresponding to the oil of these experiments) of Ma = 295 for Bo D = 0; moreover, 
the results of figure 3 indicate that the experimental results asymptotically approach 
the theoretical limit as the Bond number decreases. The transition to hydrothermal 
waves from a state of steady unicellular flow (the core of which is shown to be an 
excellent approximation to the return-flow basic state of Smith & Davis), is a very 
sharp one 

In contrast, experiments performed with ‘deep’ layers do not experience a direct 
transition from the return-flow basic state to oscillatory flow, but rather first exhibit 
a steady multicellular structure. Whether this state is the result of instability due to 
buoyancy or due to finite-geometry effects (e.g. such as the imperfect bifurcation of 
Taylor-Couette -* Taylor-vortex flow in finite annuli) is still uncertain. While even 
the smallest streamwise aspect ratio (L/d = 12) encountered in these experiments is 
comparatively large, since the present apparatus does not permit changing the depth 
of the liquid layer without also modifying the aspect ratio of the flow domain, these 
results alone are insufficient to answer this question. Recently, Cramer, Schneider 
& Schwabe (1997) have proposed a mechanism for the appearance of multicellu- 
lar structures which involves the deceleration of flow on the free surface and the 
establishment of an accompanying bulk-pressure distribution which stimulates the 
formation of multiple cells. Certainly, near the Jxrt-wall, such a rapid deceleration 
does occur in the present experiments (as shown in figured, but then one might 
question why the thinner of the two layers shown in figureX(with the higher de- 
celeration) first transitions to a hydrothermal-wave state, rather than to a steady 
multicellular state. Additional numerical and laboratory experimentation is necessary 
to completely answer this question 

Increasing the Marangoni number from a state of steady multicellular flow ul- 
timately results in the transition to time-dependent flow. In appearance, this flow 
resembles a pair of obliquely propagating hydrothermal waves riding on top of the 
steady multicellular flow. Since the primary aim of the present experiments was the 
characterization (and suppression, in Part 2) of the hydrothermal-wave instability, 
attention was restricted in the case of the oscillating multicellular state to the de- 
termination of the transitional Marangoni numbers showm in figure 3 and listed in 
table 1. 
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L Hydrothermal-wave instabilities in thermocapillary convection are known to pro- 
duce undesirable effects when they occur during the float-zone crystal-growth process, 
and perhaps in other situations. Suppression of the hydrothermal-w'ave instability 
produced in the model system of Part 1 (Riley & Neitzel 1998) is demonstrated 
through the sensing of free-surface temperature perturbations and the periodic addi- 
tion of heat at the free surface along lines parallel to the crests of the hydrothermal 
waves. 


1. Introduction 

The onset of oscillatory thermocapillary convection during the float-zone crystal- 
growth process is thought to be responsible for the appearance of undesirable stria- 
tions in material resulting from the process. This is shown convincingly in experiments 
reported by Gatos (1982), in which the heater current w'as oscillated during crystal 
growth from melt and dopant striations observed in the resulting crystal could be 
correlated with the imposed temporal oscillations. Analyses of the stability of ther- 
mocapillary convection in half-zone models of the float-zone process, such as those 
performed by Shen ef al. (1990), Neitzel el at. (1991), Wanschura et al. (1995) and 
Levenstam & Amberg (1995), yield results which identify regions of either guaranteed 
stability (energy theory) or guaranteed instability (linear theory). In the case of an 
actual crystal-growth situation, such results (were they even available) may not be 
of practical utility, since the crystal grower has little latitude in setting operating 
conditions to, say, keep the Marangoni number within a region known to be free 
of instability. In such a situation, it may be worthwhile to explore the possibility of 
suppressing oscillatory thermocapillary convection, should it occur. This provides the 
motivation for the present work. 

Riley & Neitzel (1998) (hereafter referred to as Part 1) have experimentally demon- 
strated the existence of the hydrothermal-wave instability predicted by Smith & Davis 
(1983) in thin layers of silicone oil for which buoyancy forces are quite small. This 
instability is characterized by the appearance of a single family of thermal waves 
which propagate obliquely across the layer at an angle which depends on the Prandtl 
number of the liquid. Smith (1986) has discussed the instability mechanism for this 
convective instability: for the high-Prandtl-number case appropriate to the oil used in 
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the experiments of Part 1, there is an interaction between the temperature fluctuations 
on the free surface and the bulk temperature which sustains the hydrothermal waves. 
Hence, interference with this mechanism through the modification of the free-surface 
temperature may provide a way to delay the onset of the undesirable oscillations. For 
the case of classical Rayleigh-Benard convection, Tang & Bau (1993) have demon- 
strated theoretically the possibility of delaying the onset of convection significantly 
through local modification of the temperature of the bottom wall. 

Additional motivation for undertaking an approach of this type is provided by the 
experiments of Liepmann, Brown & Nosenchuck (1982) and Liepmann & Nosenchuck 
(1982), who used periodic surface strij^neating to suppress the growth of both 
artificially induced disturbances and naturally occurring Tollmien-Schlichting (T-S) 
waves in a laminar boundary layer on a flat plate. In the latter case, the appearance 
of T-S waves was delayed significantly. These experiments were performed in water 
and the periodic heating, through the reduction of viscosity, mimicked the action of 
periodic suction at the wall. The actuation of the heating was provided by the signal 
from an upstream wall-shear-stress sensor. 

The present experiments utilize the ideas employed by Liepmann & Nosenchuck 
(1982), along with the theoretical basis provided by Smith (1986), to demonstrate 
that suppression of hydrothermal waves is possible in thin layers dominated by 
thermocapillarity. Sensing of free-surface temperature fluctuations is accomplished 
non-intrusively in real time using an infrared (IR) camera. Surface heating is delivered 
to the layer in the form of a thin sheet of infrared radiation from a COj laser which is 
out of phase with the temperature fluctuations associated with the hydrothermal waves 
at the location of the sheet. Nahas & Panton (1990) have used a similar technique to 
drive opposing flows on the surface of a jet undergoing acoustically forced capillary 
breakup, thus delaying the formation of droplets. Recently, Petrov et al. (1996) have 
employed a nonlinear control algorithm using temperature measurements near the 
free surface and heating with a thermoelectric element to suppress oscillations of 
thermocapillary convection in liquid bridges. 

2. Experimental setup and procedure 

The apparatus employed for these experiments is the same as used in the experi- 
ments reported in Part 1 ; figure 1 of that paper provides a sketch of the geometry of 
the layer. To demonstrate the feasibility of control, we concentrate on the suppression 
of hydrothermal waves for a layer d = 1 mm in depth, corresponding to axial ( L/d ) 
and transverse (W /d) aspect ratios of 30 and 50, respectively, and a dynamic Bond 
number of Bo D = 0.142. 

The output of the CO 2 laser is focused and formed into a sheet which is used 
to heat a line of small thickness on the free surface of the liquid layer downstream 
(with respect to the wave-propagation direction) of the thermocouple locations at 
which input data characterizing each hydrothermal wave are collected. The plan is 
to selectively heat the troughs of disturbance temperature, thereby locally raising the 
temperature there and hopefully disrupting the instability mechanism described by 
Smith (1986). 

The control of a particular hydrothermal wave is achieved using data measured 
specifically for that wave, necessitating the use of a delay in the feed-forward control 
scheme to be described. The signal from the infrared camera described in Part 1 
is first converted into a video signal, which is output to a video-cassette recorder 
(VCR) for recording and to a frame grabber for real-time analysis. The images from 
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Figure 1 . Overhead view of layer, showing 
(1, 2, 4 and 5) and 


hydrothermal waves and locations of data 
heating (3) lines. 



the frame grabber are processed in the following fashion, illustrated in figure 1, to 
initiate the control. Two lines (identified by the numerals 1 and 2 in figure lj) which 
are parallel to the observed hydrothermal-wave crests and located upstream of the 
known control site (line 3 in figure 1) are selected as input data locations. Along 
these lines, each of which occupies 79 pixels, the luminance data are averaged to 
smooth the noise resulting from pixel-to-pixel variations. The 79-pixel line length is 
the longest line which can be used and still perform the control in real time, given the 
limitations of the computer employed. The time histories of these averaged signals 
are used to compute the phase speed of each hydrothermal wave which passes. The 
distance between lines 1 and 2 is less than the wavelength of a hydrothermal wave so 
there is no ambiguity in the determination of this phase speed. As mentioned above, 
downstream of lines 1 and 2 is a line identified by 3 placed at the location where the 
laser heating is to be delivered. 

The output signal to the laser controller is created in the following fashion. The 
line-averaged luminance signal from line 2, called the raw input signal, is jump-limited 
and drift compensated using a running time mean. A constant input gain and a 
constant offset (which are set manually at the beginning of every experiment) are 
applied to this signal to shift it into a window (with a value range [—1, 1]) that 
is suitable for further processing and monitoring on an oscilloscope. This processed 
input signal is then temporarily stored in a buffer, so that the output signal may be 
delayed to allow for the hydrothermal-wave transit time between lines 2 and 3. Data 
read from the buffer at a location corresponding to the proper time delay are then 
inverted, gain adjusted and offset to construct an output signal within the 0-10 V 
range required by the laser controller. This output signal, tailored to each passing 
hydrothermal wave, allows the delivery of laser heating to each hydrothermal-wave 
temperature trough to suppress the oscillations. Lines 4 and 5 in figure 1 are used 
for the collection of data used to assess the effectiveness. Output gain and offset 
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may be varied manually during experiments in order to investigate the effect of the 
magnitude of laser power on the control. 

A previous attempt at control was made by Riley (1995) using a fine-wire thermo- 
couple placed near the comer of the flow field from which the hydrothermal waves 
appear to emanate. The signal coming from the thermocouple was inverted, gain ad- 
justed, phase shifted and sent to the laser controller, as above, but the phase shift of 
the signal had to be adjusted manually due to a lack of knowledge of the phase speed 
of the w'aves. The initiation of control was found to be extremely sensitive to this 
choice of phase shift. When wave suppression was achieved, it was not possible for 
it to be maintained for more than roughly a half dozen cycles, due to a combination 
of factors. First, the advection of a portion of the laser energy to the thermocouple 
location through the vigorous free-surface motion of the basic state causes a change 
in the signal there. Secondly, slight modification of the wave pattern upstream of the 
control location is observed following the successful initiation of suppression which 
necessitates a further change in phase shift of the laser input. The manual setting of 
this phase shift, coupled with the sensitivity of the system to its choice, did not permit 
this adjustment to be made, as it is automatically with the present scheme. Therefore, 
although this intrusive scheme was able to demonstrate that wave suppression is pos- 
sible for at least a few cycles, data from a single-point measurement are not sufficient 
to permit tailoring the control scheme to specific waves, as described above. 

In implementing the present control scheme described above, there are a couple 
of points to keep in mind. First, as mentioned in the previous paragraph, the free- 
surface flow is in a direction opposite to that of the streamwise component of the 
phase velocity of the hydrothermal waves and is large compared with the phase speed. 
For instance, for the case considered in this paper, the free-surface speed of roughly 
7 cm s -1 is relatively large compared with the hydrothermal-wave phase speed of 
approximately 2 mm s -1 . Hence, for the 1 cS silicone oil (Pr = 13.9) used in these 
experiments, large temperature changes created by the laser heating can be advected 
to the read lines (1 and 2 in figure 1), potentially degrading the quality of the input 
signal necessary to accomplish the suppression. Second, the hydrothermal waves which 
appear following the onset of instability are large everywhere in the flow domain (see 
also figures 8 and 9 of Part 1). On the other hand, the T-S waves of interest in the 
Liepmann &. Nosenchuck (1982) experiment are convective in nature, growing in the 
flow direction. The hydrothermal-wave instability thus appears somewhat reminiscent 
of an absolute instability in character, which may have further negative implications 
from the standpoint of suppression, particularly if the absence of oscillatory flow is 
required everywhere in the layer. 

Figure 2 show’s an overall schematic of the control experiment. The laser employed 
for these experiments is a Synrad (48-1-1 15W) CO ; laser, with variable output power 
ranging from 0 to 10 W at 10.6 jam, which is controllable with an analog input voltage 
between 0 and 10 V. The laser beam has a nominal diameter of 3.5 mm, but is allowed 
to expand in this experiment for slightly less than 2 m to assure the existence of a 
Gaussian beam. Once expanded, the beam (now 7 mm in diameter) passes through 
a 3 mm diameter aperture and is directed toward the flow apparatus by an angled 
mirror. It then passes through a zinc-selenide lens with a 250 mm focal length before 
being reflected by a concave cylindrical mirror with a focal length of -12.7 mm. This 
mirror serves to redirect the beam downward toward the free surface of the liquid 
and simultaneously expand it into a sheet. This sheet of 10.6 pm infrared radiation 
intersects the free surface of the liquid, producing the desired line-heating source for 
the control output. 



Instability of thermocapillary-buoyancy convection. Part 2 


5 



Video signal Ourput signal 


Figure 2. Schematic of the control experiments. 


For the control scheme to be successful, the thickness of the infrared sheet can 
be no greater than one-half of the wavelength of the hydrothermal waves. Since a 
1.0 mm deep layer is investigated in these control experiments, and the experimentally 
determined critical wavelength for that depth is ?. = 2.58 mm (table 2, Part 1), a sheet 
thickness of about 1 mm is desired. In order to measure this thickness, a piece of 
Plexiglas is heated and observed using the infrared camera. A graph of the line 
luminance norma! to the heated line is shown in figure 3. From this profile, the 
estimated full-width-at-half-maximum sheet thickness is 1.14 mm, which is less than 
one-half of the disturbance wavelength for the layer examined in these experiments. 
This measurement technique should yield a relatively good estimate for the laser- 
sheet thickness. In the limit of zero thermal conductivity, the image of the heated 
area yields the true sheet thickness since the heat does not diffuse away from the line 
of application. Since Plexiglas has a relatively low thermal conductivity, the thickness 
determined in this fashion is only slightly greater than the true sheet thickness of the 
infrared sheet. 

Since the C0 2 beam is invisible, the alignment of the optical system is a non-trivial 
undertaking. Optical alignment is accomplished by placing strips of Plexiglas in the 
optical path, and burning them with the beam. This allows the alignment of each 
optical component with the burned spot. Once the beam is finally expanded into 
a sheet, a rough alignment with the flow apparatus is performed using thermally 
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Figure 3. Line-luminance trace across a heated line on a Plexiglas surface. 



Wavenumbers (cm -1 ) 


Figure 4 Absorption spectrum from a smear sample of 1 cS silicone oil. 


sensitive beam-probe plates in conjunction with an ultraviolet light source. Final 
alignment of the sheet with respect to the obliquely propagating hydrothermal waves 
is accomplished using the infrared camera described in Part 1 to locate the position 
and orientation of the heated line with the hydrothermalTwave pattern present in 
the flow apparatus. Since either of two possible wave famines may appear (Smith & 
Davis^983), a mechanism allows the concave mirror to be rotated to always allow 
heating parallel to the wave crests of the observed family. 

A final note on the experimental setup and data collection concerns the absorption 
properties of 1 cS Dow Corning 200 silicone oil. Figure 4 shows an absorption 
spectrum of a smear sample of 1 cS silicone oil measured by colleagues in the Georgia 
Institute of Technology School of Chemistry and Biochemistry. Two absorption bands 
of importance are marked in the figure. The band at 3.4 pm is what is ‘seen’ by the 
infrared camera, and the line at 10.6 pm denotes the wavelength of the C0 2 laser 
corresponding to the heating. The camera is virtually insensitive to anything to the 
right of the line marked at 5.0 pm. This is important, because it allows the camera 
to detect the thermal field in the liquid, without being blinded by the much stronger 
laser emission at 10.6 pm. However, the fact that a heated line on the surface of the 
fluid layer appears as a dark stripe (see figure Sd) indicates that some of the reflected 
radiation at 10.6 pm is picked up by the IR-camera. Also note that the absorption at 
10.6 pm is relatively weak, but since this spectrum is collected for a smear sample the 
implication of this is uncertain. The data offline (1989) for 10 cS silicone oil display 
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a very similar spectrum to that measured for the 1 cS oil, and absorption results for 
a 100 pm-thick liquid sample indicate nearly total absorption. If one assumes similar 
results to hold for the 1 cS oil used in the present experiments, one would conclude 
that the heating of the liquid at the free surface actually occurs in a thin layer bsnoath at 
the free surface of roughly 100 pm thickness. For a total layer depth of 1 mm, this 
would imply that the ‘surface’ heating in these experiments occurs within the top 
10 /o of the liquid layer. We shall see in §3, however, that subsequent measurements 
indicate that this may not be the case. 


3. Results and discussion 

The control scheme described above was applied to suppress hydrothermal waves 
for a 1 mm deep layer at a laboratory Marangoni number of approximately Mai * 
575. The critical Marangoni number corresponding to this layer depth presented in 
Part 1 is Ma L = 440, so that this layer is substantially supercritical. In addition to the 
technique just described of pulsing the laser using the signal from the upstream waves, 
the influence of a constant laser power is also investigated. The effectiveness of the 
control is assessed by examining line-averaged pixel luminance data in thermographic 
images which correspond to free-surface temperature fluctuations at locations both 
upstream and downstream of the control site. For all thermographs presented in this 
paper, the cold and hot walls of the apparatus are at the left- and right-hand sides, 
respectively, and the hydrothermal waves propagate from upper-left to lower-right. 

3.1. Time-dependent heating 

This technique is the one described in §2 above, whereby the time-dependent laser 
power delivered to the IR sheet is adjusted based on the signals from the data read 
lines (lines 1 and 2 in Figure 1). The magnitude of the laser power delivered to the 
fluid layer is controlled by manually setting the values for output-signal gain and 
offset. By properly setting these two parameters it is assured that the output-signal 
troughs are very close to zero, meaning that virtually no heat is delivered to the high- 
surface-temperature locations of the hydrothermal waves. Also, by adjusting output 
gain and offset in such a fashion, the relative differences between the peak-to-peak 
and the zero-to-peak amplitudes of the output signal can be reduced to less than 
10 %. 

Figure 5 shows four instantaneous thermographs of the layer free surface taken 
while operating in this mode. The image in figure 5(a) is of the layer prior to 
initiating control, showing pure hydrothermal waves. Dark regions correspond to cold 
spots, while light regions indicate high surface temperature. In figure 5(b), control 
is attempted with the laser set at an intermediate power level. Some attenuation 
of waves is visible downstream from the heating location with hydrothermal wave 
suppression being stronger in the vicinity of the hot wall. As depicted in figure 5(c) 
this tendency is enhanced as the laser power is increased. Finally, for the power level 
used in figure 5(d), the thermograph is shown at an instant when the laser power is 
near its maximum so that the location of the heating line is also clearly visible. 

From a visual inspection of these thermographs, it is clear that the application of the 
technique has resulted in the disappearance of hydrothermal waves downstream of the 
heating location. Figure 6, however, shows evidence of the effectiveness of the control 
of a more quantitative nature. In figure 6(a), the time histories of the line-averaged 
luminance from lines 2 and 4 are plotted for a case prior to the initiation of control; in 
figure 6(b) the same quantities are shown after the start of periodic heating at a power 
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Figure 5. Infrared thermographs showing effectiveness of the time-dependent control scheme at 
various laser-aperture power (zero-to-peak) levels: (a) 0 W; (fc) 1.3 W; (c) 2.2 W; (d) 3.3 W. 


level which corresponds to the state shown in figure 5(d). It may be observed that the 
line 2 data are slightly affected (a small amplitude increase may be observed) by the 
application of periodic heating, while the peak-to-peak amplitude of the signal at line 
4 has been reduced through control from a line-averaged luminescence of roughly-fiye 
to one of 0(0.5), which is within the data noise level (determined from Fourier spectra 
to be discussed in the next paragraph). Control is achieved within roughly 30-45 s 
of turning on the laser and may be maintained virtually indefinitely As noted in 
§2, earlier experiments performed using a single-point intrusive thermocouple (Riley 
1995) to initiate the laser were only successful in obtaining a few cycles of control, 
due to a combination of a lack of knowledge of the hydrothermal-wave phase speed 
and to feedback from the heating line to the thermocouple. 

The influence of laser power level on control was investigated using an experimental 
protocol which manually stepped the laser power from zero in variable increments of 
0(0.5 V) peak-to-peak, holding the laser power at each level for roughly / minutes. 
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Figure 6 Line-averaged luminance from upstream (line 2 I dashed) and downstream (line 4 \ solid) 
data-read lines for the time-dependent control scheme: ( a ) no control (data from line 4 have been 
offset by a luminance value of 38 to allow both curves to be displayed together); (b) laser-aperture 
power (zero-to-peak) of 3.3 W, corresponding to the thermograph seen in figure 5(d) (data from line 
4 have been offset by a luminance value of 45 to allow both curves to be displayed together). 


This is more than enough time for the flow to accommodate to the new' power, 
as verified from examining the Fourier spectra computed at different time intervals 
within the thpe^-minute window (the transition time seen in the Fourier analysis is 
around 30 to 45 s). The following results are based on the amplitudes of the primary 
Fourier modes, taken during the last 37 s of these windows. Line-averaged luminance 
data are gathered at a rate of 15 Hz at each of lines 1, 2, 4 and 5 of figure 1. To 
determine the dominant component, thirty-one Fourier transforms are computed on 
a sliding data set of 256 consecutive data points. Since a discrete Fourier transform 
is employed, the observed oscillations do not consist of a pure mode, i.e. with all 
power in a single Fourier component. Rather, the vast majority of the pow-er is 
always contained in two adjacent modes which correspond to frequencies of 0.59 and 
0.64 Hz. The mean value of these two Fourier components is nearly constant during 
the 37 s period at the end of the 3-minute data window. Therefore, the analysed 
data correspond to a state, which, after the previous increase of laser power, has 
relaxed to a new equilibrium. The mean of these two largest Fourier amplitudes, is 
further time-averaged over the 37 s data-analysis window and reported as the ‘Fourier 
amplitude in what follows. Figure 7(a) shows the variation of this Fourier amplitude 
as a function of zero-to-peak voltage delivered to the laser at each of lines 1, 2, 4 
and 5 for a single experiment. Lines marked as (iHiv) on this figure correspond to 
the power levels employed in the experiments which produced the thermographs in 
figure 5 (a)-(d), respectively. At a power level slightly below 1 V, the amplitude at the 
downstream lines (4 and 5) begins to decrease, indicative of the effectiveness of the 
heating in suppressing the hydrothermal waves. 

Also observed in figure 7(a) is an upward drift in Fourier amplitude at the upstream 
lines (1 and 2). This was mentioned in the discussion of the line 2 data from figure 6. 
While the precise cause for this behaviour is not known for certain, some reasons 
may be speculated. The heating at line 3 is observed to raise the mean temperature 
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Output to IR laser (V) 

Figure 7. Spectral evaluation of the influence of laser power for the time-dependent control scheme : 
(a) single experiment, showing the behaviour at the four data-read lines (vertical lines show the power 
levels corresponding to those of figure 5); (b) line 4 data from three different sets of experiments ; 
(c) line 4 data from three different sets of experiments normalized by the corresponding line 2 data 
upstream of the heating location. 
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of the fluid at locations 1 and 2, but this has no direct effect on the values plotted 
in figure 7(a) because only the pair of dominant modes of the oscillatory component 
are being shown there (only the Fourier components corresponding to 0.59 Hz and 
0.64 Hz are considered, as discussed above). A possible indirect effect, however, may 
be important since the increase of mean temperature at lines 1 and 2 results in a slight 
increase of the time-averaged local temperature gradients upstream from the heating 
line (recall that lines 1 and 2 are near the cold wall of the apparatus). Therefore, the 
local Marangoni number Ma discussed in Part 1 is increased as laser power increases, 
and this fact could possibly explain the increases observed in the line 1 and 2 values! 
During the experiments which produced the hydrothermal waves reported here (i.e. 
in the range 450 < Mai < 570), it is found that the amplitude of hydrothermal 
waves does increase with increasing Marangoni number. A systematic investigation 
to determine the amplitude of these free-surface temperature fluctuations and their 
dependence on experimental parameters has not been performed, hence the above 
discussion remains h yp o t h eti cal . 

In figure 1(b), results of line 4 data are shown for three different sets of experi- 
ments, indicating the reproducibility of the influence of laser power level. From these 
curves, one might mistakenly conclude that the amplitude of the surface-temperature 
fluctuations at Ma L = 575 is somewhat larger than at Ma L - 570. However, this con- 
clusion cannot be drawn. The values plotted in figure 1(b) are a quantitative measure 
of the luminance variations in the analysed infrared images which depend linearly 
4l »p g nd efl o d on the surfacejtemperature variations. The linear dependence is deter- 
mined by the infrared image gain and offset, two parameters that are set manually in 
the infrared camera controller and were varied from experiment to experiment. The 
somewhat larger values plotted in figure 1(b) at Mai — 575 result from using a larger 
infrared image gain in that particular experiment. The increase of hydrothermal-wave 
amplitude with Marangoni number which was mentioned earlier therefore does not 
emerge from figure 1(b). Perhaps the most instructive and independent of assessment 
parameters is the relative attenuation of the hydrothermal waves shown in figure 7(c), 
where the line 4 data are scaled by the respective amplitudes at line 2 for each laser 
power. 

It is clear that the degree of amplitude reduction is directly correlated with power 
supplied to the laser. Further, at a power level of about 2.5-3 V, the luminance 
oscillation amplitudes have disappeared into the noise level, indicating that maximum 
possible control has been achieved. 


^•eu/«^/Ve 




3.2. Constant heating 

For this set of experimental trials, the positions of all lines in figure 1 were left 
unchanged from those of the time-dependent-heating case. The difference is that, for 
these experiments, the heating provided by the laser is constant in time, manually set 
to different values to again assess the influence of laser power. Therefore, the entire 
active-control algorithm described above is bypassed. Lines 1, 2, 4 and 5, however, 
are still used to record data for the purposes of evaluation. 

Figure 8 shows the behaviour of the dominant oscillatory Fourier component of 
the line-averaged luminance at lines 1, 2, 4 and 5 as a function of laser power for 
a single experiment. For low-to-moderate values of laser power the behaviour of 
these modes is reminiscent of what is observed in the time-dependent cases: for low 
enough values of the laser power, control is absent; as power increases, attenuation 
of the Fourier mode is seen at the locations downstream of the heating line. What 
is different in these experiments, however, is the fact that signal attenuation is also 
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Figure 8 Spectral evaluation of the influence of laser power for the constant-heating control scheme 
for a single experiment, showing the behaviour at the four data-read lines (vertical lines show the 
power levels corresponding lo those of figure 9). 


observed upstream (lines 1 and 2) of the heating line for large enough values of 
laser power, suggesting that the influence of the constant heating may be to suppress 
waves everywhere in the layer. This, however, is not actually the case, as the following 
discussion shows. 

Figure 9 shows a sequence of four thermographs of the layer corresponding to 
different magnitudes of constant laser heating; the four corresponding power levels 
are indicated on figure 8 by the letters (a), ( b ), (c) and 4 respectively. At the lowest 
of the continuous-power levels seen in figure 9(a), the heating at line 3 has no 
visible effect on the hydrothermal waves. Figure 9(h) shows a thermograph taken at a 
somewhat higher laser power level at which some suppression of wave behaviour in 
the vicinity of the heating line can be noticed. This is contrary to what is observed at 
lower power levels in the time-dependent case, where suppression is observed closer 
to the hot wall, moving toward the heating line as power is increased. Note that 
the steady output voltages employed in the constant-heating experiments are much 
smaller (by a factor of 5; see figures 7 and 8) than those (zero-to-peak) used in the 
time-dependent control scheme. As laser power is increased to the value corresponding 
to the thermograph in figure 9(c), the wave structure at locations upstream of the line 
has lost some of the regularity observed for lower power levels and a clear change 
in the propagation angle of the upstream waves is evident. Finally, the thermograph 
of figure 9(d) corresponds to a power level for which the spectral results in figure 8 
indicate a strong attenuation of waves downstream from the heating location but 
also a greatly reduced amplitude upstream. However, upstream of the heating line a 
well-defined group of waves is visible in the thermograph, but at a different angle 
than those observed for lower heating. The misalignment of these wave crests with 
read lines 1 and 2 now means that the line averages along these lines will cross wave 
crests, yielding a smaller fluctuating value, reflected in the reduced magnitude of the 
Fourier component. 

As discussed above for the case of time-dependent heating, the application of steady 
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Figure 9. Infrared thermographs showing effectiveness of the constant-heating control scheme at 
various laser-aperture power (zero-to-peak) levels: (a) 0.35 W; ( b ) 0.65 W; (c) 0.7 W; {d) 0.85 W 


heating undoubtedly affects the structure of the basic state itself. In particular, the 
mean-temperature gradient between line 3 and the hot wall is reduced by the heating; 
consequently, the local Marangoni number Ma in this region is reduced. As Ma is 
reduced, the amplitude of hydrothermal waves will decrease; if Ma falls below the 
critical value for the onset of hydrothermal waves, then waves should not appear 
there, and this is the likely reason for their observed disappearance. Therefore, it 
appears that steady power levels supplied to the layer in this manner have changed 
the underlying structure of the basic state, and hence, the instability. This is not 
the case with time-dependent heating. Evidence of this is provided by the earlier 
experiments performed using a thermocouple as the source of the input signal. As 
mentioned in §2 and §3.1, these experiments were extremely sensitive to the manually 
chosen phase shift: if this was not selected properly, hydrothermal waves covered 
the entire free surface, both upstream and downstream of the heating, in spite of 
the fact that heating was being applied. Also, following the loss of control in those 
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experiments, waves reappeared at both upstream and downstream locations. If the 
periodic heating were affecting the layer in the same manner as in the constant-heating 
experiments, one would expect the downstream waves to be absent, which was not 
the case. Consequently, the application of steady heating parallel to the crests of the 
original!) appearing hydrothermal waves is perhaps a less attractive w r ay to achieve 
control than the time-dependent scheme described in §3.1. 

3.3. Energy transmission estimate 

In addition to obtaining these results, an energy balance was performed to estimate 
the effect of the incident radiation. From the results of Pline (1989) for 10 cS silicone 
oil, it was estimated the majority of the absorbed energy at the 10.6 pm emission line 
is absorbed in a layer of roughly 100 pm thickness. Thus, the volume of fluid which 
absorbs radiation from the CO; laser is equal to the surface area which is exposed 
multiplied by the 100 pm thickness of the absorption layer. The time during which the 
heating is active is that time necessary for the local temperature depression to pass 
the heating line, which, for the case of a 1.0 mm-deep layer at the critical Marangoni 
number, is approximately Ar = 0.8 s. An approximate energy balance yields 

pc p VAT = % abs A t - mc p (T ou , - TJ, * 

L. where the term on the left-hand side represents the temperature rise of the liquid 
volume under the applied radiative flux, the first term on the right-hand side is the 
radiative energy absorbed by the layer, and the last term is the energy convected 
out of the control volume. Conductive and radiative losses from the layer have been 
neglected as they are small, and any additional convection induced by the radiative 
heating has also been neglected In the equation above, /'is the density, c p the specific 
heat, I the fixed volume beneath the applied heating (50 mm x 1.2 mm x 0.1 mm), 
Q abs is the average absorbed radiative power, m the mass flux through the control 
volume due to the basic state velocity profile, and (Tou^T in ), the average temperature 
rise of the fluid convected through the control volume. Using this result, one may 
calculate that, to effect a temperature change of AT = 0.5 C C (the peak-to-peak 
variation observed at Ma L = 575), the required peak radiative power, Q^, absorbed 
by the liquid layer should be about 25 mW. This is an extremely small value; in fact, 
it is two orders of magnitude below the applied aperture laser-power level of roughly 
3 W employed in the control experiments. 

To further investigate this apparent discrepancy, an energy-balance experiment was 
performed for a 1 cS silicone-oil layer under the same conditions experienced by 
the liquid in the flow apparatus. An insulated liquid layer of volume V = 10 cm 3 
was placed at the location normally occupied by the flow apparatus and exposed to 
the CO; laser sheet at a power level comparable to that employed for the control 
experiments (P * 2.5 W @ the laser aperture) for a period of 5 minutes. The change 
in the mean temperature of the layer following this exposure was AT » 0.9 °C. 
Accounting for the optical losses in the experimental setup, it was estimated that the 
amount of energy reaching the free surface of the liquid is roughly 15% of the 2.5 W 
output by the laser. Therefore, of the roughly 1 10 J of energy incident upon the 
liquid layer, only 12.5 J were actually absorbed, indicating an absorption efficiency 
of 1 1 %. Using the calculated absorption efficiency of 1 1 %, the aperture laser power 
required to deliver the estimated peak absorbed radiative power of Q ptak se 25 mW is 
approximately P , 1.5 W, which compares reasonably well with the observed values 

of Pi * 3.0 W used in the control experiments. 

It is interesting to note that the measured absorption efficiency of 11% for the 1 cS 
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silicone oil is in contrast to the assumption made on the basis of results for 10 cS 
silicone oil, which Pline (1989) determined to be totally absorbing at a layer depth 
of 100 pm. Therefore, the smear sample of the 1 cS silicone oil shown in figure 4 
probably provides a correct picture when it indicates that the 1 cS silicone oil used in 
these experiments is not very absorptive at the 10.6 pm spectral line corresponding 
to the laser output. 

From the experimental results, one observes that excellent control of the hydrother- 
mal waves is achieved with an output voltage to the laser of 3 V zero-to-peak (see, 
e.g. figure 7c). Using the efficiency calculated above, the power absorbed in the fluid 
layer for this output voltage is roughly 50 mW. This flux is approximately ten times 
that calculated as necessary to erase the hydrothermal-wave free-surface temperature 
depression. 


4. Conclusions 

The experiments described above have demonstrated the feasibility of suppressing 
oscillatory thermocapillary convection in thin liquid layers subject to a horizontal 
temperature gradient. Based upon the instability mechanism suggested by Smith 
(1986), control was achieved by heating troughs of low disturbance temperature with 
a sheet of IR radiation. The active feed-forward control scheme employed in these 
experiments is conceptually simple but effective for the suppression of hydrothermal 
waves in a high-Prandtl-number liquid. 

Two features of the described experiments were key to their success: (i) the non- 
intrusive sensing of the surface-temperature fluctuations using an infrared camera; 
*-nd (ii) the real-time determination of the hydrothermal-wave phase speed, allowing 
the control input to be precisely tailored to each hydrothermal wave. Earlier attempts 
at control by Riley (1995) using a thermocouple inserted below the free surface of 
the liquid achieved limited success due to the reliance on a single-point measurement 
and, perhaps, to the intrusive nature of the thermocouple inserted through the liquid 
free surface. The manual adjustment of the phase and amplitude of the forcing in 
those experiments made it impossible to suppress the oscillations for more than a few 
cycles. 

There remain questions regarding the effectiveness of the C0 2 laser in providing 
surface heating to the 1 cS silicone oil employed for these experiments. Clearly, the 
laser is sufficiently powerful to accomplish the task, but the precise depth range of 
the applied heating is still unknown, as is the importance of this range. The work 
of Smith (1986) implies that, for high-Prandtl-number liquids, surface-temperature 
disturbances are coupled, through a convective process, to thermal disturbances in 
the bulk to sustain the instability. Hence, heating supplied to the bulk as well as to 
the surface may also provide a means of breaking this coupling. 

Finally, the eventual use of such a technique in suppressing oscillations occurring 
during the float-zone crystal-growth process requires additional thought. First, such 
crystal-growth situations involve liquid metals of extremely small (0(1O -2 )) Prandtl 
number, so that the mechanism of the instability is different, as observed by Smith 
(1986) for planar layers and later by Levenstam & Amberg (1995) for a cylindrical 
geometry. In the planar case, hydrothermal waves propagate transversely across 
the layer, normal to the mean temperature gradient, meaning that the elimination 
of oscillations at the cold wall (solidification front) would require the simultaneous 
suppression of waves in the entire layer. For a cylindrical geometry exhibiting relatively 
few waves propagating in the azimuthal direction, this may or may not be possible. 
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